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We propose a configuration of a magnetic microtrap which 
can be used as an interferometer for three-dimensionally 
trapped atoms. The interferometer is realized via a dynamic 
splitting potential that transforms from a single well into two 
separate wells and back. The ports of the interferometer are 
neighboring vibrational states in the single well potential. We 
present a one- dimensional model of this interferometer and 
compute the probability of unwanted vibrational excitations 
for a realistic magnetic potential. We optimize the speed of 
the splitting process in order suppress these excitations and 
conclude that such interferometer device should be feasible 
with currently available microtrap technique. 
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I. INTRODUCTION 

Since the first realization of magnetic traps [0J|] and 
guides with current-carrying conductors on a chip, 
a large variety of magnetic potentials have become ex- 
perimentally accessible, which would be impractical or 
even impossible to realize with macroscopic coils. The 
splitting of two-dimensionally trapped atom clouds has 
been demonstrated J5|||, and recently, we were able to 
split and unite a three-dimensionally trapped cloud of 
rubidium atoms in a chip trap 0. 

Current experiments aim at populating single quan- 
tum states of such microtrap potentials with either an 
atomic ensemble (i.e., creating a Bose-Einstein conden- 
sate), or indeed with a single atom. One promising ap- 
plication of such a system would be an integrated atom 
interferometer on a chip ||. The small size and mono- 
lithic construction of such a device suggests its suitabil- 
ity for "real-word" applications. Moreover, the fact that 
magnetic potentials may be "engineered" on the chip en- 
ables novel interferometer schemes with features quite 
different from more traditional atom interferometers [...]. 
Here we study a scheme in which the particle wave of a 
single, trapped atom is coherently split up and reunited 
by a time- varying magnetic potential (fig. 1). Splitting 
occurs in one dimension, while tight confinement in the 
remaining two dimensions leads to an effective ID situ- 
ation. This is in contrast to || where the dynamics of 
the splitting is in two dimensions. As depicted in fig. 1, 
interference occurs between the lowest two vibrational 
states, \ifo) and of the splitting potential (the inter- 
nal atomic state remains unchanged). A phase-changing 
interaction in one "arm" (stage H in fig. 1) translates into 
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FIG. 1. Scheme of the trapped atom interferometer: one or 
several atoms are prepared in the vibrational ground state of 
the single well potential (I) . When the well separates, the wave 
function evolves adiabatically into a symmetric delocalized 
state (II a). The phase of the wave function in each potential 
well can be assumed to evolve independently and monitors 
sensitively external electric and magnetic field gradients (II 
b). As the potential wells reunite, the antisymmetric state 
transforms into the first exited vibrational state whereas the 
symmetric one retransforms into the ground state (III). 



a change of the relative populations in \ipo) and \ipi) when 
the potential is recombined. As in other interferometers, 
a longer duration of stage II leads to a larger accumu- 
lated phase (i.e., a larger arm length). However, unlike 
the situation in most free-atom schemes and the guided- 
atom scheme proposed in [||, in our scheme the atom 
does not move nor does its wave function spread during 
this stage: the propagation along the traditional interfe- 
rometer path is replaced by the evolution in a constant 
potential (stage II), which leaves the position and the 
physical size of the wave function unchanged. This in- 
terferometer is thus particularly well suited to measure 
local fields and interactions, which presents an advantage 
over experiments with propagating atomsf]. One could, 
for instance, measure the phase shift arising from a two- 
body collision |l0| or the amount of decoherence induced 
from a nearby surface p| . 

In this paper, we present a detailed analysis of this in- 



1 There is a subtle difference between atom interferometers 
with beams and with trapped atoms: in spatial beam splitters 
atoms are slowed down when the energy of the transverse 
state increases. There are currently studies on the way how 
this effect can be explored for enhanced detection schemes of 
the outgoing state 
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terferometer scheme, employing a realisitic magnetic po- 
tential which can be implemented with currently avail- 
able microtrap technique. We consider the case of an 
individual trapped atom, a situation which is also tar- 
geted by experiments under way (for a study of a BEC 
in an idealized ID potential, see jL2j). The potential 
is created by the simple conductor configuration shown 
in fig. 2. The current Iq together with the external field 
.Bo, y provides tight confinement in the yz-plane. The two 
currents 7 ex t together with the homogeneous field compo- 
nent B x close this 2D trap in axial direction, completing 
the single- well potential (fig.||a, left). The current I c cre- 
ates an adjustable "bump" in the center of this trap, and 
thus induces the splitting. Increasing I c transforms the 
potential from single- well to double- well (fig.||a, right), 
in loose analogy with the first passage through the beam- 
splitter of a Michelson-Moreley interferometer for light. 

To achieve a good fringe contrast, it is essential that 
no higher-lying vibrational states be excited during this 
splitting. Therefore, the crucial part of the interfe- 
rometer is the quantum dynamics during the splitting 
and merging process. The splitting (merging) of the 
wave functions occurs as the quantum states adiabati- 
cally evolve in the varying potential. We analyze these 
dynamics in a one-dimensional model, using analytic ex- 
pressions for the microtrap magnetic field. We numeri- 
cally determine the energy eigenstates of 87 Rb atoms in 
the given potential and then trace the dynamics of an 
initial state, using the eigenstates as a time-dependent 
basis. 

We show that successive vibrational levels in the initial 
trap evolve into pairs of degenerate states when the po- 
tential is split. In this "sensing state" , the wave functions 
are composed of two identical oscillator states in the left 
and in the right well. Either of the two parts can acquire 
a phase shift independent of the other one, reflecting e.g. 
an additional small field gradient or the presence of an 
additional atom in one of the wells. When the poten- 
tial is transformed back into a single well, the population 
of the vibrational levels depends on the phase difference 
that is picked up in the degenerate states. Fig.|T] illus- 
trates this process: first, the system (i.e. one or several 
atoms) is prepared in the vibrational ground state. Upon 
separation of the potential, this state evolves into a sym- 
metric state that spreads over the two potential wells. In 
an analogous manner, the antisymmetric first vibrational 
level transforms into an antisymmetric delocalized state. 
As the system's Hamilton operator is symmetric through- 
out the whole process, it cannot induce transitions be- 
tween states of opposite symmetry, and the eigenstates 
can always be chosen of well-defined parity. 

If the symmetric and antisymmetric state are spa- 
tially separated far enough, they degenerate, and the left 
(right) localized state can be constructed as sum (dif- 
ference) of the symmetric and antisymmetric state. A 
perturbation of the potential, which does not have even 
parity, will lift this degeneracy in favour of the localized 
states. These localized states make up for the classical in- 
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FIG. 2. Layout of the interferometer conductor pattern. 



terferometer arms, measuring very sensitively deviations 
from an ideal symmetric potential or interactions with 
other atoms. 

In the following section, we investigate the separation 
process using the ID-potential taken from the microtrap 
device sketched in fig. |[ We establish the quantum me- 
chanical equation of motion and use first order pertur- 
bation theory to determine the amount of vibrational 
excitations. Assuming a linear variation of the current 
I c , we find the excitation probability lower than 2% if 
the separation takes 60 ms or longer. This indicates that 
an experimental realization should be possible, and the 
situation can still be improved when an arbitrary varia- 
tion of I c is allowed. We therefore dedicate section III to 



a method which minimizes non-adiabatic excitations by 
finding the most appropriate time dependence for the 
shape of the potential (here controlled via I c ). Such 
method is of interest not only for interferometers. It ap- 
plies to all cases of time-dependent potentials, and it can 
even be transferred to spatially varying potentials such 
as beam-splitters. For our interferometer, this method 
helps to reduce the splitting time by a factor of two, 
at the same time reducing the excitation probability by 
more than a factor ten. 



II. THE TRAPPED ATOM INTERFEROMETER 

The microtrap device that we propose for the interfe- 
rometer is a symmetric arrangement of wires as depicted 
in fig.^j. Its potential is similar to the one that we used 
in the merging experiment with thermal atoms Q], but 
it is scaled down to a wire distance of 20 fim and sim- 
plified to produce a strictly symmetric potential. The 
quantum state computations are made for 87 Rb atoms in 
the \F = 2, nip = 2) ground state, the effective potential 
being U[x) « h ■ 1.4MHz • B(z)/G. 

The current Iq in the central wire and the homoge- 
neous field -Bo, j/ create a two dimensional quadrupole field 
which strongly confines the atoms in the j/z-plane. Each 
of the crossing wires contributes a longitudinal field mod- 
ulation of Lorentzian shape (see p3|): The two currents 
^oxt together with the field component B^ x generate two 
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FIG. 3. a) Shape of the magnetic splitting potential for 
characteristic values of the control parameter s (see eq. (^) 
and (§) ), b) eigenstates of 87 Rb atoms in the specified poten- 
tial. 
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FIG. 4. Energy eigenvalues of the system's Hamilton oper- 
ator as the trapping potential is divided into two wells. Neigh- 
boring states of opposite symmetry form pairs and degenerate 
as the potential wells separate. 
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valleys along the longitudinal axis, which do not appear 
separate if the trap is located far enough from the sur- 
face. The current I c with its direction opposite to the 
two external currents is used to split the Ioffe-Pritchard 
potential into two neighboring wells (fig.|]a). Choosing 
the parameters as 
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the trap is located 35 ^m above the surface, yielding a 
transversal oscillation frequency of Wtrans ~ 27r 53.7kHz. 
The parameter s determines the shape of the trap, run- 
ning from for one single well to 1 for separated wells. 
The point s — 1 has been chosen such that the two low- 
est vibrational levels of each well are clearly separated 
(i.e., the two lowest sets of states are both degenerate). 
The time-dependence of the system's Hamiltonian is ex- 
pressed via the function s(t). In a simple approach, s 
may be chosen to vary linearly in time, but as we will 
discuss in section II], an optimized function s(t) can be 
found which minimizes vibrational excitations during the 
splitting (merging) process. 

In fig.Ha, the resulting magnetic field along the longi- 
tudinal axis e x is displayed for characteristic values of s; 
the transverse potential minimum is plotted against the 
longitudinal position. The plots below show the eigen- 
states of 87 Rb atoms (\F = 2, nip = 2)) in this field as 
they are numerically computed from the Schrocdinger 
equation (eq. (j^) below). The corresponding energy 
eigenvalues, measured relative to the minimum value of 
the potential, are given in fig.||. 

For s = 0, the four lowest levels correspond to the 
states of a harmonic oscillator with quantum number 
n = 0, . . . , 3 and oscillation frequency w s= o ~ 27r 190 Hz. 
We use the quantum number n to identify the eigenstates 
as \ip n (s)) throughout the whole evolution. As the value 
of s is raised, the vibrational levels evolve into symmet- 
ric and antisymmetric delocalized states. At s = 1, the 
four lowest levels form two sets of degenerate states, their 
energy being E 2 k = E 2 k+i = huj s=1 (k + |), k = {0, 1}, 
u> s= i w 2-7T 240 Hz. At each stage, the separation of the 
transverse levels (> 50 kHz) is much larger than the sep- 
aration of longitudinal states involved. For this reason, 
the longitudinal states do not intermingle with the trans- 
verse levels, even if the system's symmetry is slightly dis- 
turbed. The quantum dynamics is therefore adequately 
described by a one-dimensional model. 

In order to make the interferometer work properly, the 
atomic wave function should follow ideally the (time- 
dependent) eigenstates \ipk(t)) of the system. If the po- 
tential is varied too fast, the evolution is non-adiabatic, 
i.e. vibrational excitations are generated. For the inves- 
tigation of these excitations, we will focus on the first 
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half of the interferometer cycle: we use a time-dependent 
interaction picture to compute the time scale on which 
the separation process can be lead adiabatically. 

The (time-dependent) basis for the computation is 
found by solving the time-independent Schrodinger equa- 
tion 



H(s)\ip k (s)) = huj k (s) \<fk(s)) 
with the Hamilton operator 



H(a) 



P 

2m 



Mb 9f nip |B(s,r)| 



(6) 



(7) 



where s takes the role of a mere parameter. For the given 
magnetic field, the eigenfunctions have been computed 
numerically and are displayed in fig. ||b. 

The natural phase evolution of the eigenstates can be 
included into the basis and yields the ansatz 



m)) 



(8) 



The equation of motion for the coefficients Ck(t) is ob- 
tained when eq. (||) is inserted in the time-dependent 
Schrodinger-equation with the Hamiltonian (Q)q 



■{<p k (t)\-\ vn {t)) 



(9) 



Given that a single eigenstate \(fi(t = 0)} is prepared in 
the beginning, and further assuming that the transition 
probability into other vibrational states is small, first or- 
der perturbation theory can be used to determine the 
coefficients c/ and the corresponding transition proba- 
bilities Pif. 



Cf(t) 



S «r(«/M— MK^ctjii^ct))* (io) 



Pif(t) = \c f (t)\ 



(11) 



The coupling (<p f {t)\&\tpi(t)) = <¥>/(*)! £l¥><(*)> ft to 
higher levels is directly proportional to the rate f at 
which the control parameter s is changed. Therefore, if 
all levels are separated by a minimum energy htoo, the 
transition amplitudes can be made negligible by choos- 
ing an appropriate duration for the process. Conversely, 
if at certain instants some energy levels degenerate, this 
will create large transition amplitudes unless the coupling 
coefficient (<Pf(s)\-^\<fi(s)) between these levels vanishes 
at the points of degeneration. In the trapped atom inter- 
ferometer presented here, we encounter such degenerate 



2 The time-dependence of ujk(i) and \(fik(i)} is explicit 
through the control parameter s: Wfc(t) = Wfc(s(t)) etc. 
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levels. But as the states that degenerate are of opposite 
symmetry throughout the complete evolution, the cou- 
pling coefficient remains zero for all times. Therefore, 
the excitation probability can be made arbitrarily small 
by choosing the process duration long enough. 

This consideration is confirmed by numerically evalu- 
ating expressions ( |l0| ) and ( |TT| ) for either of the interfe- 
rometer levels (\ipo),\ipi)). In a first approach, the sepa- 
ration parameter has been chosen linear in time s = t/T. 
Fig. H shows the transition probabilities into the neigh- 
boring interferometer levels which contribute largest to 
all vibrational excitations. The data indicate that the 
excitation probability is less than 1% if the separation 
process takes longer than 60 ms. 

This is an encouraging result, as it seems experimen- 
tally realizable. Moreover, as the time-dependence of the 
potential can be freely chosen, one can adjust the speed 
of the separation process f in order to further reduce 
excitations. This is of general interest, because a linear 
variation of the control parameter s s not necessarily the 
best choice. Indeed, one wishes to find a method that 
optimizes the process irrespective of its parametrization. 



III. OPTIMIZATION 

In this section, we develop a scheme which minimizes 
vibrational excitations in time-dependent potentials. In 
a slight variant, this method can equally be used to find 
an adequate shape for a beam-splitter potential. 

In order to optimize the adiabaticity of the separation 
process we first take a look at the coupling term from eq. 

(<Pf(t)\f t \Vi(t)) = (M^\Vi(s)) J , (12) 

which is proportional to the process speed 4f and to the 
coupling coefficient a(s) = (<Pf(s)\j^\<Pi(s)). 
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Intuitively, one can increase the process speed if 
a(s) is small, and decrease it in the opposite case. Fur- 
thermore, the process speed should be adapted to the 
energy difference of the levels involved, being the more 
increased the further the energy levels lay apart from each 
other. Last not least, one has to avoid discontinuities in 
the process speed including the start and the end of the 
separation. In the following, these intuitive rules will be 
substantiated into a set of differential equations to yield 
an optimized process control s(t). 

We assume that the process is lead during < t < T 
and that the separation parameter at t — T is s(t — T) = 
1. Indeed, we want to fix a shape of the control param- 
eter s which does not depend on the process duration. 
Therefore, we implicitly assume that s(t) = s(t, T) can 
be written as 



fl (t,T) = *(~ 1) 



(13) 



The goal is then to fix some maximum excitation proba- 
bility e 2 and to find an appropriate shape for the function 
s(t,T) which minimizes T a di a b fulfilling the condition 



\c f (T)\ <e VT>T adiab . 



(14) 



If, by some chance, the distance of energy levels 
Aw (s(t')) = ujf — uji is constant throughout the process, 
the transition amplitude ct(T) appears as the Fourier 
transform of a(s) 4f : 



ds 







t a{s)-dt 



(15) 

for Auj = const. (16) 



If, in addition, a(s) happens to be constant over the pro- 
cess, the solution of the problem is simple: the shape of 
the process speed 4| should be chosen such that it pro- 
duces the least amount of side bands possible in a Fourier 
transformation. An appropriate shape would e.g. be a 
Blackman pulse p4[ : 



ds 
di 



1 
T 



— cos(2tt-) + — cos(4tt-) 



(17) 



which can be directly integrated to yield s(t). 

The idea of the Fourier transform can be extended to 
the more general case. A substitution of the time variable 
t by some new variable r can be made in a way that the 
argument of the exponential in eq. (|l^) becomes linear 



J\u{t')dt'=^rr{t) 



(18) 



and that r runs from to 1 during the process. The time 
scale Tq will be part of the optimization result. Equation 



@ then assumes the form of a Fourier transform of 
some new expression u(t) [J 



r 1 t 

Cf (T) = / e^ T u{r)dT 
Jo 



ds T 

u(t) = a{s)-^- 



(19) 
(20) 



The expression u(t) is a generalized coupling term, 
acting in the transformed time frame r. As above, one 
can now choose a shape for this coupling term u(t) (how- 
ever, not its amplitude) and will obtain the probability 
amplitude as its Fourier transform. 

After the optimization strategy is chosen, it remains to 
solve the equations (|l8| ) and (|20|). One might be tempted 
to deduce the relation J£ from eq. ( p"8| ) and insert it into 
eq. ( po|) to solve directly for s(t). Unfortunately, this re- 
sults in an intractable problem. Instead, one can take ad- 
vantage of the substitution already made and first solve 
for s t (t). The relation between r and t is then estab- 
lished in a second step. This way, the problem is split 
into two differential equations the first of which gives the 
amplitude of u(t), and the second of which determines 
the time scale To used in the substitution. These two val- 
ues determine size and scale of the probability amplitude 
c f {T). 

The first differential equation involves the shape of the 
function that is chosen for the generalized coupling term 
u(t), and it is a direct consequence of eq. (|2C 



ds T u{t) 
dr a(s T (r)) 



(21) 



It is important to note, that although the time is 
used to shape the coupling term it(r), its amplitude does 
not correspond to the overall process speed. Instead, 
the amplitude of u(t) has to be adjusted such that the 
solution matches the boundary conditions s(r = Q) = 0, 
s(r=l) = 1. This can for instance be done by iteratively 
solving eq. ( |2l| ) for different amplitudes of u(t). 

The second differential equation establishes the rela- 
tion between t and r and arises from the substitution of 
t (eq. ( |T§| ) ), once that s T (r) has been determined: 



(22) 



Choosing T = To, this equation can be solved numeri- 
cally, and one finds To as the point in time, for which 
r(t) reaches its boundary t(Tq) = 1. 

The result for the transition amplitude is now com- 
pletely described by equation ([l9|), the amplitude of u(t) 



3 In the following equations, the index r marks the fact that 
the functional dependence of the parameter s is on r, not on 
t. 
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and the time scale resulting from the choice of the pulse 
shape. The optimized evolution of the control parameter 
is computed from the concatenation of s T (r) and r(t): 



s(t) 



(23) 



If this optimization is applied to the trapped atom 
interferometer, the probability for non-adiabatic excita- 
tions can be considerably reduced. Fig. || shows the ex- 
citation probabilities for a process speed which has been 
optimized to suppress the transition |<^>o} — * With 
the optimized control, the separation can be done within 
30 ms, thus reducing the complete interferometer cycle to 
60 ms with an overall excitation probability of less then 
10~ 3 . 

These parameters suggest that an experimental real- 
ization of the scheme is indeed feasible. It remains, of 
course, a difficult task to prepare the atoms in the ground 
state and to detect the atoms selectively in different vi- 
brational states. However, this seems achievable using 
Bose-Einstein condensates of low density. Another is- 
sue is stability against gradients of magnetic stray fields. 
In our case, the sensing states of the interferometer lie 
~ 6 /mi apart. During a sensing time of 60 ms, a gradi- 
ent b x = would lead to an additional dephasing of 
A$ « 2-7T 50 • b x /—. A suppression of stray gradients to 
less than 1 mG / cm would therefore reduce the dephasing 
to A(/> < §§. 



find parameters which suggests an experimental imple- 
mentation in the near future. Grounding on the theoreti- 
cal results, we have developed an optimization scheme for 
the reduction of vibrational excitations that is indepen- 
dent of the system's parametrization. Applying the opti- 
mization to our interferometer potential, we have found 
a cycle of duration T — 60 ms with excitation probability 
less than 10~ 3 . 
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IV. CONCLUSION 



In conclusion, we have studied a dynamic potential in- 
terferometer working with three-dimensionally trapped 
atoms. We have used a time-dependent interaction pic- 
ture to describe the quantum state evolution and we have 
computed probabilities for non-adiabatic transitions into 
neighboring levels. For a realistic magnetic microtrap we 
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